Enhancement of gaps in thin graphitic films for 
heterostructure formation 



J. p. Hague 

The Open University, Walton Hall, Milton Keynes, MK7 6AA, United Kingdom 
Abstract. 

There are a large number of atomically thin films with similar structure to 
graphene. These films have a spread of bandgaps relating to their ionicity, and 
also to the substrate on which they are grown. Such films could have a range 
of applications in digital electronics where graphene is difficult to use. I use the 
dynamical cluster approximation to show how electron-phonon coupling between 
film and substrate can enhance these gaps in a way that depends on the range 
and strength of the coupling, as well as the length scales of spatial fluctuations. 
One of the driving factors in this effect is the proximity to a charge density wave 
instability, which could in principle generate gaps in graphene. The enhancement 
at intermediate coupling is sufficiently large that spatially varying substrates 
and superstrates could be used to create heterostructures in thin graphitic films 
with position dependent electron-phonon coupling and gaps, leading to advanced 
electronic components. 



PACS numbers: 7f.45.Lr, 71.38.-k, 73.22.Pr, 73.61.Ey 

1. Introduction 

The 2D material graphene has made headhnes over the past decade for its remarkable 
properties. Often overlooked is the availability of other two-dimensional graphitic 
materials which have a similar structure and properties to graphene, but with a 
direct bandgap that is lacking in suspended graphene. These gapped compounds have 
the potential to make graphene compatible digital transistors, semiconductor lasers 
and solar cells. Atomically thick graphitic materials with honeycomb lattices and an 
inherent direct bandgap formed because of strong ionicity include boron nitride (BN) 
[T] (band gap 5.6eV [2J) and other materials can be grown in very similar hexagonal 
Wurtzite layers, such as InN (band gap 0.7-0.8eV) InSb (0.2eV [1], possibly down 
to 45meV on certain substrates f5]), GaN ( 2.15eV [6]), and AIN (6.28eV ^7]). It 
has also been reported that small gaps can be formed with a similar mechanism in 
graphene-gold-ruthenium systems [5] and graphene-SiC systems (there is some debate 
about the latter piTU]). 

Recently, I used a self-consistent mean-field theory to show that gaps in atomically 
thin hexagonal materials may be modified by introducing strong electron-phonon 
coupling through a highly polarisable superstrate [11] [12]. Similar interactions 
between graphitic monolayers and substrates form polaronic states and affect the 
overall electronic structure of the monolayers, as shown by quantum Monte Carlo 
simulations for highly doped thin graphitic films [13]. Effective electron- electron 
interactions can be induced via a strong coupling between the electrons in atomically 
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thick monolayer and phonons in a highly polarisable substrate because of limited out of 
plane screening, similar to that seen for quasi-2D materials [Hj . Such interactions have 
been experimentally demonstrated between carbon nanotubes and a Si02 substrate 
|15| . and are necessary to account for the lower mobilities of graphene on Si02 |16| . 

Here, I make calculations beyond the mean-field theory by using the dynamical 
cluster approximation (DCA) to compute the effects of electron-phonon interaction 
on electrons in atomically thick graphitic materials, where a gap has been opened 
because of ionicity. I present results computed with a high order perturbation theory 
consistent with Migdal's theorem (which allows neglect of vertex corrections for low 
phonon frequency and weak coupling) and discuss the effect of long range interactions 
and spatial fluctuations. 

Besides the use of electron-phonon interactions with substrates, the possibility 
of tuneable gaps has mainly focused on graphene. Following a theoretical proposal 
[TTl [T5] . bilayer graphene has been observed to have a gap that can be tuned 
by applying an external electric field [ini HOj- Electron confinement in graphene 
nanoribbons leads to gaps |21| . and high quality nanoribbons can be made by 
unzipping nanotubes [22] or using patterned SiC steps [23]. Very wide bandgaps 
can be formed by functionalizing graphene with hydrogen (graphane) [211 [HJ [55] and 
fluorine (fluorographene) [571 [55]. 

This paper is organised as follows: A model Hamiltonian for the interactions 
between graphitic monolayers and substrates is introduced in Sec. [2] The perturbative 
expansion and dynamical cluster formalism used to solve the model are discussed in 
Sec. [3] Sec. |4] presents details of gap enhancements and the spontaneous formation 
of a charge density wave state. A summary and conclusions are presented in Sec. [5] 

2. Model Hamiltonian 

The Hamiltonian required to describe the motion of electrons in thin films with 
honeycomb lattices has a basis of two atoms. Typically, electron motion within the 
plane is described using a tight binding model, and ionicity is taken into account 
with the potential ±A on the two sublattices. With a highly polarisable substrate, 
there is additional electron-phonon interaction between the electrons in the film and 
phonons in the substrate, which may be long range (i.e. momentum dependent). A 
Hamiltonian with these properties has the form, 

H = Hth + Hci-ph + Hph (1) 

where Hth is the tight binding Hamiltonian representing the kinetic energy of the 
electrons in the monolayer, Hd-ph describes the electron-phonon interaction, and 
Hph is the energy of the phonons in the substrate (treated as harmonic oscillators, 
and including both kinetic and potential energy of the ions). 
The tight binding part of the Hamiltonian is written, 

+ A{al^aka ~ cl^Cka))- (2) 
The first part represents the kinetic energy, where 0^ = — i exp(ifc.^j) and 
6i are the nearest neighbour vectors from A to B sub-lattices, <5i — a(l,-\/3)/2, 
S2 = a(l,— •\/3)/2 and S3 = (— a,0) and a is the spacing between carbon atoms 
in the plane (the tilde is used to avoid confusion with the creation and annihilation 
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operators). Electrons with momentum k are created on A sites with the operator a\ 
and B sites with c\^. The second part represents the interaction between electrons in 
the monolayer and a static potential, either induced by the substrate (in the case of 
graphene) or by ionicity (in monolayers of III-V semiconductors). Here, A sites have 
a higher potential, A and B sites are lower in energy by -—A. Breaking the symmetry 
between A and B sites in the bipartite honeycomb lattice gives rise to a gap. 
The phonon part of the Hamiltonian is, 

H^i, = Y,%{b% + d\d^) (3) 
q 

where phonons with momentum q are created on A and B sites with &^ and 
respectively. Typically, the phonon dispersion, f2q is taken to be momentum 
independent as a good approximation to optical phonons. Typical phonon frequencies 
vary from 10s to 100s of meV. For example, in BN phonon energies range from llOmeV 
for transverse acoustic phonons at the K point to 200meV for optical phonons [29]. 
Due to ionicity, sites have a net charge, so strong coupling between electrons and 
phonons is expected. 

Finally, the interaction between electrons in the monolayer and phonons in the 
substrate (or superstrate in the case of graphene on a substrate) is, 

kq 

+ II 5q [aLq«fc(4 + d-q) + 4_qCfc(fet + b-q)\ (4) 
hq 

where the coupling constants g and g represent interactions between electrons and 
phonons on the same sub-lattice and different sub-lattices respectively. A lattice 
Frohlich electron-phonon interaction with the form, 

gmin) oc exp(-|n - m\/Rsc)[{n - rnf + c^]"^/^, (5) 

has been proposed for layered quasi-2D systems [H]. Experiment has demonstrated 
that this form explains interactions between electrons in carbon nanotubes placed 
on Si02 [E], and is necessary to account for reduced mobilities when graphene is 
placed on a substrate [16]. The screening radius, Rsc controls the length scale of the 
interaction, c is the distance between the graphitic thin film and surface atoms in 
the substrate. In the following, I take (? = 2a^, since the distance between graphene 
and substrate (which are typically bound by van der Walls interactions) is likely 
to be slightly larger than between the very strongly bound carbon atoms in the 
graphene layer. Ionic, graphitic materials may bind more strongly to appropriate 
ionic substrates leading to shorter c, which in this work is represented (in combination 
with screening effects) with a reduced Rsc- In practice, the effects of changing c and 
Rsc on the form of the effective electron-electron coupling mediated by phonons are 
very similar. 

The model used here has some similarities to the ionic Hubbard model [30]. In 
the ionic Hubbard model, the ionicity (introduced by an analogous parameter A) acts 
against the Mott insulating state (which is caused by the repulsive Hubbard U). In 
contrast, in the model here, the parameter A acts with the electron-phonon coupling 
to form a charge density wave (CDW) insulating (gapped) state. 
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Figure 1. (Colour online) Dynamical cluster approximation (DCA) sub-zones for 
cluster sizes of up to Nc = 100. Within each sub-zone, the self-energy is taken to 
be momentum independent. This allows Green functions to be calculated in the 
thermodynamic limit, and convergence properties are particularly good as Nq 
is increased. Note that the only symmetry taken into account is translation in 
fc-space. 



3. Method 

The electron-phonon Hamiltonian described above is extremely difficult to solve 
exactly using numerical methods. An approximate solution can be made using 
perturbation theory within the dynamical cluster approximation formalism. The 
dynamical cluster approximation (DCA) [211 is one of the possible ways of 
extending dynamical field theory [|] (DMFT) [55] so that it can be applied accurately to 
low dimensional systems. The Mermin-Wagner-Hohenberg theorem indicates that the 
significant non-local fiuctuations found in one- and two- dimensions are expected to 
eliminate spontaneous long range order, potentially rendering results from mean-field 
theories qualitatively wrong [331 131]. DCA resolves this problem by re-introducing 
spatial fluctuations. 

When applying DCA, the Brillouin zone is divided up into Nc sub-zones centred 
about a momentum vector Ki (see Fig. [I]) consistent with the symmetry of the 
whole system. Within each sub-zone, the self-energy is approximated as a momentum- 
independent function, so the Green function can be coarse grained by integrating over 
the sub-zone, 

G{Kui)= [I(z + //) + A^3-*fc-S(fs:i,z)]-i. (6) 



X Where DMFT is used to approximate low dimensional systems, it is often known as the local 
approximation 
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In finite size techniques, the number of particles is related to the number of 
momentum points used in the calculation of the self energy. In contrast, the 
DCA coarse-graining step involves an infinite number of momentum points, so the 
thermodynamic limit is satisfied for any cluster size. In the context of the perturbation 
theory for the Migdal-Eliashberg theory used here, DCA has particularly good 
convergence properties in cluster size Nc, so in principle smaller clusters can be used 
leading to a significant improvement in computational efficiency |36|, I37j . 

In several previous studies, Nc = 4 and Nq = 16 DCA clusters have been used 
to understand Hubbard interactions on hexagonal/triangular lattices (see e.g. [55]'). 
I briefly discuss the subzone schemes for hexagonal lattices with larger Nc- For the 
lattices used here, the simplest way of defining the sub-zone vectors is: fei = kx/\/Nc, 
k2 = k2/y/Nc with the vectors Ki = nki + mk2 with n and m integers. Here, the 
reciprocal lattice vectors are fci — (27r/3a, 27r/-\/3a) and fc2 — (27r/3a, — 27r/A/3a). 
There are likely to be other valid lattices that can also be used, where the lattice and 
sub-lattice are oriented at different angles. However, the lattices used there are the 
simplest to implement. 

Even for the simple cases considered here, the resultant lattices can be ordered 
into groups. Clusters with Nc = (3n)2 {Nc = 9, Nc = 36, Nc = 81 etc.) have 
sub-zones centred on the K and K' points (here n > 1 is an integer). Those with 
Nc = (3n — 1)^ {Nc = 4, Nc = 25 and N^. = 64) make a second set where 3 sub- 
zones share a corner at the K and K' points, and the third set of Nc = (3n -I- 1)^ 
{Nc = 16, 64, 100 etc.) where 3 sub-zones share a corner at the K and K' points and 
an edge with the full Brillouin zone. Since the self energies would be identical in the 
3 zones around the K and K' points in the latter 2 cases, they will poorly describe the 
physics at the K and K' points (which is especially important for graphene). This is 
why I use only the (3n)^ series. 

To establish which k point belongs to a sub zone, it is sufficient to find the closest 
Ki point corresponding to the centre of the sub-zone (subject to shifts of reciprocal 
lattice vectors). The edges of the shapes defined in this way are the hexagons in the 
figure. 

The perturbation theory used here can be seen in Fig. [2] Panel (a) shows the 
Hartree diagram. For the symmetry broken states, this cannot be absorbed into the 
chemical potential, and is the main contributor to modification of the gap. The Fock 
potential (shown in panel (b)) is responsible for frequency dependence of the self- 
energy. Phonon propagators are modified using a Dyson equation (panel (c)) which 
modifies the phonon frequency, and can lead to further enhancement of the band gap. 
Therefore, 
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Figure 2. Perturbation theory used in the work presented here, (a) Hartrce 
diagram. For the symmetry broken states, this cannot be absorbed into the 
chemical potential, and is the main contributor to modification of the gap (b) 
The Fock potential is responsible for frequency dependence of the self-energy 
(c) Dyson equation for the phonon propagator. This renormalizes the phonon 
frequency, which can lead to further enhancement of the band gap. 



where Q represent the centres of the coarse-grained cells for phonon momenta. 

D;j\Q) = d-^{Q)-Ui,iQ) (11) 

where 

T 



ILij{Q,uj,) = -2-— 2^ G^-jiQ + K,ujn +LOs)GjiiK,ujn) (12) 

and the non-intcTacting phonon propagator is, 

d,,{Q,iOs) = XijiQ)ny{fl^ + ujj). (13) 

Here the dimensionless, momentum dependent electron-phonon coupling is 
defined using the following procedure. First, the coupling for a single k value is 
defined to be A^'^) = |/f )|V2Mt02 ^^ere fl^^^ = fi"""^ = Eie''=-«^/fl„ ft^ 
j2^^ik-R,+iaf^^^.^^ and /f = E,'5*'"'''~''/R.-ia. Here, /r, = go{R,)^2Mn/h, 
the lattice vectors are Ri = hai + fha^, a\ = 3i/2 + •\/3j/2 and a2 = 3i/2 — •\/3j/2. 

The dimensionless electron-phonon coupling can then be related to \j{Q) via, 

\ (ij) 

4. Results 



X,,{K,)=2NcX ^'''^^ ^ (14) 



The aim of the work presented here is to use dynamical cluster approximation (DC A) 
to examine how charge density wave (CDW) gaps in graphitic thin films vary with 
electron-phonon coupling. I take T = Q.02t, A = O.lt and Q = O.Olt. Noting 
that t is typically on the order of an eV, these values correspond approximately to 
room tempcaature, phonon frequencies, fl, of 10s of meV and A a few hundred meV, 
consistent with thin films of materials such as InSb, or reported gaps in some graphene 
on substrate systems. 

I start by computing self-energies to show the relative contributions of on- and off- 
site terms, and the effects of varying the interaction range. The computed self energies, 
including the relative contributions of Hartree and Fock diagrams for A = 2 and 
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Figure 3. (color online) Comparison of Hartree and Fock parts of self energy, 
S = S' + iS". The main contribution to the self energy is the real part of the 
Hartree term, which is momentum independent. N.B. The Fock term is still the 
most important contribution in some cases, for example the inverse mass depends 
on derivatives of the self energy, so this will be given by the Fock term. T = 0.02t, 
A = O.lt, Q = O.Olt, A = 2 and A^c = 9- 



Nc — 9 are sho-wn in Fig. [3j The real part of the Hartree diagrara is momentum and 
energy independent. It is the largest part of the self energy, and directly contributes 
to the enhancement to the gap, and also is the main contributor to spontaneous CDW 
order. N.B. The Fock term is still the most important contribution for some properties, 
for example the inverse mass (not considered here) depends on derivatives of the self 
energy, so this property -will be given by the Fock term. The real parts of the on-site 
self energies are significantly smaller for the long-range Frohlich interaction than for 
the Holstein model. Other differences are that the off-site self energy is small for 
the Holstein interaction (N.B. it is not zero because of effective off-site interactions 
mediated through the phonon self-energy), "whereas it is of similar magnitude to (but 
smaller than) the onsite self-energy for the Frohlich interaction. 

The main aim of this paper is to understand the role of electron-phonon coupling 
range in the enhancement of gaps. Figure |4] sho-ws ho^w the enhancement varies "with 
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Figure 4. (color online) (a) Gap enhancement for a Holstein interaction 
with A/t = 0.1, comparing results from the dynamical mean-field theory 
and DCA. There are only small corrections due to momentum dependence, 
(b) Gap enhancement for the long range Frohlich interaction, Rsc = 1 (c) 
Gap enhancement for the long range Frohlich interaction, Rsc = 3 (d) Gap 
enhancement for the long range Frohlich interaction. The Nq = 144 point at 
A = 2 in panel (d) demonstrates that convergence has been achieved. Here 
T = 0.02t and f2 = O.Oli. 




Figure 5. (color online) Spontaneous symmetry breaking. Here T = 0.02t, 
n = O.Olt, N(j = 1 and A = 0. Note that a larger A value is required to see CDW 
for longer range interactions, but the effect is still present. 



interaction range and with cluster size. Larger cluster size means longer range spatial 
fluctuations with length scale ^ V^c- Figure |4|^a) shows gap enhancement for a 
Holstein interaction with A/t = 0.1, comparing results from the dynamical mean- 
field theory and DCA. N.B. The enhancements will be smaller for systems with large 
ionicity [TTJ[T2]. There are very small corrections due to momentum dependence. As 
the screening radius, R^c increases, the enhancement decreases. Spatial fluctuations 
have no effect on this set of diagrams for the Holstein interaction. Figs. Qb) and 
[4j[c) show results when — 1 and Rg^ = 3 respectively. Essentially, the effective 
A decreases with R^c- This is most extreme for Rsc oo and Fig. |4|^d) shows 
gap enhancement for the long range Frohlich interaction. The reductions in gap 
enhancement due to spatial fluctuations are also most pronounced for long range 
coupling. Convergence is found for Nc ~ 81 indicating that spatial fluctuations extend 
over approximately 9 sites. Clearly the key to finding a strong effect is to reduce the 
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effective range of the interaction, or alternatively to increase A. This may be achieved 
by choosing substrates that sit closer to the graphene, or increasing screening from 
the substrate. 

The phonon self-energy plays an important role in the gap enhancement. Fig 
|4ja) also shows the enhancement effect when Nc = 1 if the polarisation bubble (PB) 
is neglected and the curve can be compared with the full theory with Nc = 1. The 
phonon self energy augments the enhancement by increasing the value of the phonon 
propagator at the Brillouin zone centre. Thus, increasing the electron self-energy, 
which is proportional to the phonon propagator. 

It is also interesting to determine if gaps can spontaneously form from the 
electron-phonon interaction. This is explored in Fig. [5j which shows spontaneous 
charge density wave (CDW) symmetry breaking when Nc = 1. For sufficient A, a 
CDW state can be found for all values of Rgc (this is not visible in the figure for large 
values of Rsc because A values are too small). While it is expected that the CDW 
state will become more difficult to form as cluster size increases (owing to the Mermin- 
Wagner-Hohenberg theorem [331 EHj) the result shows that the full system with A 7^ 
is on the cusp of a CDW state, and this is why the gap is strongly sensitive to the 
electron-phonon coupling. 

5. Summary and conclusions 

In summary, I have investigated gap formation and enhancement in a model of 
atomically thin graphitic materials. Electron-phonon coupling and range has been 
varied, and the effects of higher order corrections to the phonon propagator have 
been considered. The effect of reintroducing fluctuations around the mean field limit 
has also been investigated using the dynamical cluster approximation. Higher order 
corrections to the perturbation theory increase the gap enhancement. It is found that 
gaps are enhanced by electron-phonon interactions for all interaction ranges, with the 
enhancement decreasing as interaction range increases. The re-introduction of short 
range fluctuations is also found to cause further decreases to the enhancement. In 
order to see these effects in real materials, it would be necessary to reduce the range 
of the interaction, either by choosing ionic substrates that sit close to the graphitic 
monolayer and have a strong coupling with electrons in the fllm, or by promoting a 
strong screening. 

One of the driving factors of this enhancement is the proximity to a charge 
density wave state for a material without ionicity (A = 0) such as graphene. I 
have shown that sufficiently large coupling between electrons and phonons can lead 
to spontaneous CDW order in the mean-field solution. This instability to order 
shows why there are significant gap enhancements at large coupling when ionicity is 
introduced. While a spontaneous CDW state would be unstable in a pure 2D material 
because of the Mermin-Wagner-Hohenberg theorem [Ml [3S], a quasi-2D material 
with weak hopping between layers should be able to generate a CDW state at the 
temperature characterised by the interplane hopping. Only small interplane hopping 
due to tunnelling is necessary to stabilise the CDW state (to stabilise a state below 
T ~ 500if, inter plane hopping on the order of 50meV, or around 1% of the in plane 
hopping is required). Owing to the spontaneous symmetry breaking, an appropriately 
layered heterostructure of graphene and a wide gap insulating material such as BN 
might be expected to generate small spontaneous gaps of useful size due to CDW 
formation. 
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Figure 6. Schematic showing the possible use of two substrates with different 
electron-phonon coupUng to thin graphitic films to make a heterostructure in 
a single thin film of a graphitic material. This could be manufactured by 
laying down an interface between the two substrate materials before cleaving 
perpendicular to the interface and then depositing the thin film. 

The results here suggest that an interesting possibility would be to use the 
electron-phonon interaction to make position dependent changes to the bandstructure 
of the thin film (for example by adding a spatially dependent superstrate with phonons 
that strongly couple to electrons in the thin film), a method that is potentially easier 
to control than trying to deposit neighbouring thin films with interfaces in the plane. 
The proportional gap enhancements predicted here arc similar to those between GaAs 
and AlGaAs [31], so it is plausible that thin film hctcrostructurcs or quantum dots 
could be built up in this way (see Fig. [6|. Another possibility would be to tune 
inherent gaps in III-V semiconductors with the electron-phonon interaction, so that 
they become optimal for applications such as solar cells where the efficiency is highly 
sensitive to the gap size. Clearly graphitic thin films warrant further study to assess 
their full capability for novel electronics. 
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